A scalable, tunable qubit, based on a clean DND or grain boundary D-D junction 
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Q\ ■ Unique properties of a ballistic DND or grain boundary D-D junction, including doubly de- 

t-H \ generate ground state with tunable potential barrier between the "up" and "down" states and 

» ! . non-quantized spontaneous magnetic flux, make it a good candidate for a solid state qubit. The role 

of quantum "spin" variable is played by the sign of equilibrium superconducting phase difference 
^5 , on the junction, which is revealed in the direction of spontaneous supercurrent flow in equilibrium. 

Possibilities of design-specific simultaneous operations with several integrated qubits are discussed. 
The pronounced shift from "software" to "hardware" in theoretical research on quantum computing (QC) |l] [1]] 
is a good measure of growing confidence that QC can be realized on practically interesting scale (of at least lCr 
qubits). Though first experimental realizations of QC used such technologies as NMR and ion trapping, the problem 
of scalability for such approaches still looks formidable. Therefore much effort is directed at search for a practical 
solid state qubit (SSQ), with natural candidates being such mesoscopic devices as quantum dots , mesoscopic 
^ \ Josephson junctions and superconducting single-electron transistors (parity switches) H . The evident advantage 
$H ■ of a SSQ is scalability, where all the potential of existing solid state technologies could be used while among the 
problems the main are (1) to achieve quantum beatings between distinguishable states of a single qubit, (2) to prevent 
loss of coherence during calculations, and (3) to minimize statistical dispersion of the properties of individual qubits. 

The problems (1) and (2) pose specific difficulties in an SSQ due to huge number of degrees of freedom coupled to 
it, and to necessity to fine-tune two states of a mesoscopic system chosen as working ones to a resonance. 

A possibility to circumvent these obstacles is presented by Josephson systems with d-wave cuprates, which violate 
time-reversal symmetry and as a result have doubly degenerate groundstate |^,^| with a potential for quantum beatings 
(macroscopic quantum tunneling), or at least quantum noise M. Ioffe et al. recently incorporated this property in 
their "quiet qubit" design which uses tunneling (SID) or dirty SND junctions with equilibrium phase difference 
§ : ¥>o = ±§. 

i—i" Let us consider a device shown in Fig. la. Its main part is a clean mesoscopic D-D junction (i.e. ballistic DND or 
D-(grain boundary)-D junction). Though we will concentrate on the DND qubit design, the same considerations are 
applicable, mutatis mutandis, to high quality grain-boundary D-D junctions, where recently an analogous current- 
phase dependence was observed [||. The grain boundary region where the superconducting gap is suppressed, is 
naturally modeled by a normal conductor with the same lattice parameters and chemical potential as in the super- 
y—{ [ conducting banks, which only enhances the amplitude of purely Andreev scattering in the system. 

■ The terminal B of the junction is formed by a massive d-wave superconductor; in a multiple-qubits system, they 
will all use it as a common "bus" bar. The terminal A is small enough to allow - when isolated - quantum phase 
fluctuations. It is essentially the sign of the superconducting phase difference tp between the terminals A and B that 
plays the role of " spin variable" of quantum computing. The collapse of the wave function is achieved by connecting the 
terminal A with the external source of electrons ("ground"), thus blocking the phase fluctuations due to phase- number 
uncertainty relation ||. So far, the best way to do this is presented by using a "parity key", PK (superconducting 
single-electron transistor) , which only passes Cooper pairs, and only at a certain gate voltage V g . Other parity 

keys, with different parameters, are used to link adjacent qubits, allowing for controllable entanglement (Fig. lb). 
Such an architecture allows reasonably easy way to integrate a large number of qubits in a ID or 2D matrix. We will 
see that it also provides a natural way of preparing all qubits on the same bus in the same initial quantum state, thus 
J-^ ' facilitating the implementation of error correction algorithms. 

The readout of the state of a qubit is simplified by the presence of small spontaneous, non-dissipative currents 
and spontaneous fluxes (of order 10 -2 — 10 -3< I>o depending on the setup) concentrated in the central part of the 
DND junction which have opposite directions in two degenerate equilibrium states. While too small to lead to 
unwanted inductive coupling between the qubits or decoherence, they can be still used to read out the state of the 
qubit once it was collapsed in one of the states with ±ipo, e.g. using the magnetic force microscope tip M (which is 
removed during the computations). Collapsing and reading processes are thus time separated, and the computation 
results will be automatically preserved for the time limited only by the thermal fluctuations. 



> 

o 



ON 
-i— > 

a 

I 



- 1—^ 

X 



* email: zagoskin@physics.ubc.ca 



1 



A necessary condition for a qubit to work is t t ^ unneling) < t g(ate app u caU on) < td(ecoherence)- Here a clean DND 
junction has a very important advantage following from the fact that the absolute value of equilibrium phase difference, 
| ipo\, can now vary from to tt depending on the angle between the crystal axes of the d-wave superconductors and 
the ND boundary fll2| . In practice, since the lattice structure of the d-wave superconductors allows only a limited 
set of easy cleavage directions, e.g. (010) and (100), the equilibrium phase can be varied by preparing a steplike ND 
interface, with the equilibrium phase determined by the relative weight of Andreev zero- and tt— levels (coupled to 
the lobes of d-wave order parameter in A and B with the same or opposite sign respectively) , produced by such an 
arrangement. Since the shape of the effective potential barrier between states with ±<^>o depends on ipo, it can now be 
varied. Therefore the tunneling rate can be chosen in exponentially wide limits to achieve the optimal performance. 
(In an SND junction the equilibrium phase can also be varied, but it cannot be made less than ir(y/2 — l)/\/2) jl^.) 
Besides, this allows to fix the working interval of the device to \ip\ < tt, since due to exponential dependence of the 
tunneling amplitude on the barrier action, for |<^o| < tt/2 the tunneling to the states in the next cell, ipo — > ipo ± 2tt, 
can be completely neglected, whatever the inductance of the system. 

The supercurrent through the normal part of the system is carried by a set of Andreev levels formed by reflections at 
the ND boundaries |l3|]l4f | . We will calculate it using the quasiclassical approach following from Eilenberger equations 
p5[ which is well suited to our problem. Here we can consider "Andreev tubes" along the quasiparticle trajectories in 
the normal part of the system []l2"| , [l6 17 each carrying the supercurrent density which in our case should be written 
as 



' It sinhp£[r, n\/l T 



P =i 



(1) 



Here the critical current density j c ~ {evp n)/(Aj? S), S and n being area and half-perimeter of the system. Equation 
(|I|) is valid in the limit £o "C S/H, the characteristic separation between the superconducting electrodes and is thus 
applicable to the case of superconducting cuprates. (The opposite limit of large coherence length is considered in 
|18|.) The quantity C[r, n] is the length of the quasiclassical trajectory linking two superconductors which passes 
through the point r in the direction n; </ 3 [£[ r J n ]] is the phase gain along this trajectory (in this paper we do not 
concern with the effects of current-induced field, and therefore in the absence of external fields this is simply the 
phase difference between the superconductors on the ends of the trajectory, including the extra tt if C[r, n] connects 
the lobe of the d-wave order parameter with opposite signs jt2)). The normal metal coherence length It = v F /2iTkBT , 
and limp takes into account effects of weak elastic scattering by nonmagnetic impurities (in the ballistic regime, by 
definition, li mv 3> S/Tl) . We have used standard approximation of steplike behaviour of the order parameter at the 
ND boundary, and neglected the own magnetic field of the supercurrent |l9| . 
The total supercurrent density at a point r is thus given by 



j(r) 
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Calculating the total current flowing in A, we find in the limit 
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1 + Z(Q) 
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where F((p) is the 27r-periodic sawtooth of unit amplitude, and the imbalance factor Z(Q) determines the equilibrium 
phase difference Q 



\fo\ 



1 - Z(0) 



(4) 



In the setup of Fig. la 
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The current-phase dependence (||) and corresponding Josephson energy Ej((p) — ^ J v dipl(ip) are plotted in Fig. 2, 
and current density distribution in the normal part of the system is shown in Fig. 3. The vortex pattern is clearly 
seen. The spontaneous flux in the system is 



$ s ~ k(0) 



j c II 2 K(Q)eN ± v F 1 k(Q) v f 

C C 137 TT c 
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where k < 1 is a geometry-dependent attenuation factor (e.g. in SND junctions k = by symmetry if ND boundary 
is parallel to (100) or (010) |^@); Nj_ ~ H/Xf is the number of transport modes in the system. 

Let us make some estimates. Taking the size of the system ~ 10 3 A, vp ~ 10 7 cm/s, we find that $ s ~ k ■ 10 _3 $o- 
The magnetic moment of the spontaneous current will be of order m s ~ K(N±eHvF)/c ~ k-IO 5 //^. The tunneling rate 
between the states is T ~ ujq exp[— U(0)/?uvo], where the frequency of oscillations near ±<po is ujq ~ yJN±_eeQ/h and 
the height of the potential barrier U(0) oc (<pQ/2n) 2 m. Due to spatial Andreev quantization, there are no elementary 
excitations in the normal part with the system with energies below e ~ Hvf/211 ~ 10 -15 erg. At temperatures below 
T = e/ks 10K thermal excitations are frozen out, and dissipation can only be due to interlevel transitions generated 
by ac Josephson voltage generated by phase fluctuations. Therefore it will be absent if 



2e < Vj hyj < Cp 2 > ~ huo < e, (7) 

which can be rewritten as €q < e/N± (where the charging energy cq = 2e 2 /C, and C is the capacitance of the 
terminal A), or ujq < Vp/211. The latter condition is a physically clear requirement that the quantum oscillations 
of superconducting phase allow time for readjustment of Andreev levels in the system (which is indeed ~ 2IL/vp, 
the time necessary for the electron and Andreev reflected hole to travel across the system). Otherwise the coherent 
transport through the normal part of the system cannot be established, and dissipative currents flow instead. The 
maximum value of ujq allowed by the above limitation is w max ~ 10 s . (That is, the capacitance of the terminal A 
cannot be lower than C m i n = 2e 2 N±/l ~ 10 _11 F.) The corresponding tunneling rate is 

r max ~ w mex exp[-N±{ip /2%) 2 ]. (8) 

Therefore we would require ipg ~ 0.2ir to have tunneling rate in the 100 MHz region. 

For a system of integrated DND qubits of Fig. lb the bulk d-wave "bus" provides, the possibility for operations 
performed over all the qubits simultaneously by creating a supercurrent flow along the bus. In particular, one can 
easily prepare the whole register of qubits in the same (up or down) state. This is an attractive property, e.g. for 
implementation of a quantum correction algorithm. If take the size of a unit qubit with its periphery as ~ 5 • 10 3 A, a 
2D 100 x 100-qubit block will occupy only an area of 50x50 /xm 2 , which is realistic to keep below the dephasing length 
due to thermal excitations. 

Application of quantum gates to individual qubits can be effected in various ways, lifting the degeneracy between 
up/down states either by directly applying localized magnetic field to a qubit (using a magnetic scanning tip), by 
creating local supercurrents in the bus, or by using laser beams with circular polarization. The entanglement of the 
states of adjacent qubits is achieved simply by opening a key between them for a certain time. 

In conclusion, we have suggested a new design for a solid state superconducting, scalable qubit. Besides using the 
degeneracy of the ground state, common to all D-D junctions, it strongly relies on unique properties of clean DND 
or grain boundary junctions: tunability of the equilibrium phase difference across the junction, and spontaneous 
currents and fluxes in equilibrium. The former allows to optimize the design in order to achieve the fastest possible 
tunneling rate (which is vital in order to beat the dephasing processes) and thus makes it easier to integrate qubits 
in a computer. The latter presents an easier way to manipulate and read out the state of a qubit. Our estimates 
show that there is a real chance to create a working solid state qubit along these lines using existing experimental 
possibilities. 
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FIG. 1. (a) Superconducting DND qubit: A,B are d-wave superconductors, N normal conductor, PK parity key, M scanning 
tip, £1 the mismatch angle between the lattices of A and B. The cut in B is here along (110) and (llO) directions. Positive 
lobes of d-wave order parameter are shaded, (b) Multiqubit register. Terminal B plays the role of the bus bar. (c) Version of 
(b) using grain boundary (G) junctions. 
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